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1 Introduction 

Chaos is a typical property of many-dimensional nonlinear systems. Its role 
it revealed in various problems of astrophysics and cosmology. Chaos made 
to revise the two-hundred year old views on the evolution of Solar system. 
Theory of interstellar matter, dynamics of star clusters and galaxies at present 
cannot be considered without chaotic effects. 

Astronomical topics themselves had remarkable impact on the develop- 
ment of chaotic dynamics. The Hcnon-Hciles system, one of the first systems 
with revealed chaotic properties, was proposed for the study the motion of 
a star in a galactic potential. Much earlier, Poincare's classical work on the 
foundations of the theory of dynamical systems emerged from the problem of 
small perturbations in the planetary dynamics. 

In the present lectures I will discuss only several astrophysical and cos- 
mological problems. The choice of the problems is determined with the aim, 
first, to cover as broad topics as possible, second, to show the diversity of 
approaches and mathematical tools. I will start from planetary dynamics, 
moving to galactic dynamics, to cosmology, and to the instability in the 
Wheeler-DeWitt superspace. For pedagogical reasons I will describe the tech- 
niques, such as the estimation Kolmogorov-Sinai entropy, the dealing with 
hyperbolicity in pseudo-Riemannian spaces, so that they can applied for any 
other problems. Obviously, numerous other problems, methods and results 
remain out of these lec tures, m ost of them, however, can be trac ed from the 
references; for chaos seJ^EEl f or applications of our interest se JSEEEl 

I will start from a brief review of the elements of theory of dynamical 
systems, to introduce the main concepts used in the subsequent chapters. 
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2 Elements of Ergodic Theory 



2.1 Dynamical systems 

Ergodic theory is the metric theory of dynamical systems, i.e. which deals 
with spaces for which a measure is defined but not a metric. 

In the following brief account of elements of smooth ergodic theory, we 
will concentrate on the classification of dynamical systems by the degree of 
their statistical properties; for details se tr l u l I 

The key concept is obviously, that of the dynamical system. Initially 
the dynamical systems were understood as mechanical systems, however later 
that term was generalized to variety of physical systems of non-mechanical 
origin. Cosmological solutions of Einstein equations can be considered as 
such examples. 

Dynamical system (M, B, fx, T) is considered defined if M is a smooth 
manifold, B is a cr-algebra of measurable sets on M, and fi is a complete 
measure on B, and T* is a one-parameter group of diffeomorphisms defined 
by vector field v 

One-parameter groups are called flows, by the term borrowed from hydro- 
dynamics, and below we will give elements of the classification of flows. The 
apparent abstractness of the definition implies quite general and natural prop- 
erties for physical systems. 

2.2 Classification of dynamical systems, mixing, relaxation 
Flows are called ergodic, if for any measurable invariant set A 

T l A = A = T-*A, (2) 
its measure y, takes only the values 

M(A) = |J. (3) 

One can show that for measure-preserving ergodic flows the time-average 
almost everywhere equals the phase space average 

/ fdn= lim - [ f(T- T (x))dT. (4) 
Jm t^ 00 t Jo 

In physical literature this property is often considered as a definition of an 

ergodic system since it is enough and sufficient. The property of ergodicity is 
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one of rare definitions of smooth ergodic theory which can be generalized also 
for spaces with infinite measure. Ergodicity, however, is a weak statistical 
property and therefore is less important for actual physical problems. 

The far more importance for statistical physics of another property, mix- 
ing, has been established firstly by Gibbs. 

Ergodic theory provides definitions for mixing of various degrees. 

Weak mixing is indicated by the condition for V/, g € L 2 



lim — 



f( T T x)gd^i- I fdfi / gd\x 

M JM JM 



dr = 0. (5) 



The 'weakness' of the property of weak mixing can be seen from the 
following limit 

1 I* 

lim- / | fi(T~ T A f)B) - n(A)fj,(B) \dr = 0, (6) 

t^oo t J Q 

implying that T l A becomes independent of the set B only if some parts of 
the trajectory are not taken into account. 

Note the absence of the factor T/t' and hence increase in the convergence 
rate in the definition of the property of mixing 



lim / f(T t x)gdn= [ fdfi f gd^. (7) 

VM JM JM 



Analogically the property of m-fold mixing for m functions is generalized as 
follows 

lim / f f 1 (T^x)...f m (T^+ -+ t ^x)d l i = TT / fidfi. (8) 

fi,..,t m ^ooj M r=o"' M 

These properties describe systems with increasing statistical properties in the 
sense that, systems with mixing possess the property of weak mixing, and 
those with n-fold mixing also that of mixing and weak mixing but not vice 
versa. 

Systems with mixing are evidently also ergodic ones. However, for the 
systems with mixing, as opposed to ergodic ones, a set A £ B evolves in such 
a way (preserving its measure and connection) that the measure of the part 
which intersects the set B € B tends in time to be proportional to the measure 
of B 
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Compare this limit with the following one for ergodic systems 

1 f* 

lim- / n(T- T An B)dr = n(A)n(B), (10) 
t^oo t J 

The latter limit is said to converge in the Cesaro sense, while the limit for the 
mixing case is ordinary converging. In other words, for ergodic systems the 
initial fluctuations tend to zero only in the time-average, for mixing systems 
their absolute value decreases as well. 

Hence the property of mixing guarantees the existence of a final state of 
measure fi to which 

ln{A)=iio(T t A) 1 AeB (11) 

tends smoothly, so that 

lim f fdnt = f fdfi. (12) 

In physical terminology the final state is called \i equilibrium, the process of 
tending to that state is the relaxation. 

Even stronger statistical property, K-mixing, possess K-systems (Kol- 
mogorov systems) systems for which the following limit 

lim sup | fj,(A (IB) - n(A)n(B) \= (13) 

t — >oo 

exists, where the upper limit is taken for the smallest er-algebra containing 
T to A for t < to. Kolmogorov systems possess n-fold mixing of arbitrary n. 

Strongest statistical properties are possessed by hyperbolic (Axiom- A), 
Anosov, and Bernoulli systems. 

We will define Anosov system^^l w hich by their statistical properties are 
equivalent (isomorphic) to Bernoulli shifts. 

The flow /* is of Anosov type if for its all trajectories {/*} there exist 
subspaces E s (f t (x)), E u (f t (x)) of the tangential space TMft^, and numbers 
C > 0, A > 0, such that 

TM fHx) =E s (f(x))®E u (f(x)), 

drE s (f{x)) = E s (f t+T (x)), drE u {f(x)) = E u (f t+T (x)), 

and for alH > 0, one has 

|| dfv \\< Ce~ xt \\v ||, veE s ; 
II dfv \\> C-V* || v ||, veE u . 

a As mentions Smale^^lhe had found the horseshoe transformation, the classic example of 
Axiom-A systems in Brazil, on Rio beaches. 
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The subspaces E s and E u are stable (converging) and unstable (expanding) 
subspaces. 

For physical systems this definition implies exponential instability at each 
point of the phase trajectory and at any small perturbation. Anosov systems 
are subclass of hyperbolic systems and they possess an important property 
of structural stability. Roughly it means that the perturbed systems possess 
the property of the unperturbed one, i.e. the strong instability acts towards 
preserving of the properties of the system. Though strictly speaking the 
conditions of Anosov systems are never or almost never are satisfied for real 
physical systems, nevertheless it appears that the structural stability can be 
peculiar to certain types of strongly instable physical systems. 

Geodesic flow on a compact manifold with negative constant curvature 
is an example of an Anosov system, and was studied long ago by Hadamard, 
Hopf and Hedlund. Their works had inspired 

KryknEI 

to apply those ideas 

to physical systems. 

If the systems with mixing can tend to equilibrium by any law (e.g. poly- 
nomial), hyperbolic systems tend to that state exponentially. 

2.3 Kolmogorov-Sinai entropy 

The problem of distinguishing different features of dynamical systems, and the 
formulation of corresponding characterizing criteria is a central one in ergodic 
theory. Much efforts in this direction were concentrated on the study of spec- 
tral properties of dynamical systems, until in 1958, Kolmogorov discovered 
the new metric invariant, the entropy. 

Consider the entropy of a splitting & of the measurable manifold M 



d 




(14) 



i=l 



where (iGS and 




(15) 



d 




(16) 



i=i 



Then the Kolmogorov-Sinai (KS) entropy h is the limit 

fc(/) = sup lim -H(C), 



(17) 
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where 



71-1 



c = V 



(18) 



3=0 



and the upper limit is taken over all measurable splittings. 

Dynamical systems with positive KS-entropy h > are usually called 
chaotic, while those with h = are called regular ones. In particular, Anosov 
and Kolmogorov systems, which are typical systems with mixing, have positive 
KS-entropy h > 0, while most of only ergodic ones have h = 0. Therefore the 
latter are not considered to be chaotic according to this definition. 

For the above mentioned geodesic flows on spaces with constant negative 
curvature R < 0, the KS-entropy equals 



The KS-entropy is related to the Lyapunov characteristic exponents Aj 
via the Pesin formula 



we see that a system with at least one non-zero Lyapunov exponent has pos- 
itive KS-entropy. The use of Lyapunov exponents for many-dimensional sys- 
tems is not always well defined, nevertheless it was efficiently applied for stellar 
systemic 

Finally let us mention another important characteristic of dynamical sys- 
tems, the correlation function, defined by 



Although at present estimates of the correlation functions (including numer- 
ical results on some billiards) exist only for a few dynamical systems, for 
Anosov systems it has been shown that the correlation functions decay expo- 




(19) 




(20) 




(21) 



nentially, i.e., 3a gg ',f3,t > so that 



K g '(t)\ <a g>gl exp(-/3i), 



(22) 



where 



/3~/i(/). 



(23) 
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3 Chaotic Solar System 



otitic 



effects in planetary dynamics. For detailed reviews I would refer to' 
where various evidences of chaos, particularly in the asteroid belt, in the 
motion of comets, are discussed, along with the methods of overlapping res- 
onances and estimation of Lyapunov exponents, Wisdom-Holman symplectic 
mapping and other techniques used in those studies. 

Before considering the stability of the Solar system, let us formulate the 
two key theorems, Poincare's and Kolmogorov's, which were crucial in the 
efforts on this long-standing problem. 

N-dimensional system is considered as integrable if its first integrals 
I\,...,In in involution are known, i.e. their Poisson brackets are zero. As 
follows from the Liouville theorem if the set of levels 

M I = {I j (x) = I?,j = l,...,N} 

is compact and connected, then it is diffeomorphic to N-dimensional torus 

T N = {(9 1 ,...,9 N ),modd2Tr}, 

and the Hamiltonian system performs a conditional-periodic motion on Mi. 
Poincare theorem states that for a system with perturbed Hamiltonian 

H(I,<p,e) = H {I) + eHi(I,<p,e), (24) 

where I, cp are action-angle coordinates, at small e > no other integral ex- 
ists besides the one of energy H = const, if Ho fulfills the nondegeneracy 
condition, 

det\du/dl\ ^ 0, (25) 

i.e. the functional independence of the frequencies to — dHo/dl of the torus 
over which the conditional-periodic winding is performed. 

Though this theorem does not specify the behavior of the trajectories of 
the system on the energy hypersurface, up to 1950s it was widely believed 
that such perturbed systems have to be chaotic. 

Kolmogorov's theorenff^ f 1954 ; the main theorem of Kolmogorov- 
Arnold-Moser theory, showed that at certain conditions the perturbed Hamil- 
tonian systems can remain stable. 

It states: 

If the system (24) satisfies the nondegeneracy condition (25) and H\ is 
an analytic function, then at enough small e > most of non-resonant tori, 
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i.e. tori with rationally independent frequencies satisfying the condition 



^ nu k ^ 0, 



(26) 



do not disappear and the measure of the complement of their union set 
fx(M) -> at e -> 0. 

KAM-theory was initially considered as supporting the views on the sta- 
bility of the Solar system, though it says nothing about the limiting value of 
the perturbation e. 

However, though the level of direct applicability of the KAM theory for 
the Solar system remains not clear, it appears that, the joint application both 
of theoretical and numerical methods at present computer's possibilities is 
rather efficient. 

The frequency map technique developed by Laskar | 20 | 21 | 22 |23| j g ^ asec j 
on the approach of KAM theory. This method enabled numerical treatment 
of long-term planetary evolution in terms of a perturbed Hamiltonian system 
using the idea that, if a quasi-periodic function is given numerically on the 
complex domain, then it is possible to approximate it via a quasi-periodic 
function with an accuracy higher than that given by standard Fourier series. 
Namely, the quasi-periodic fun ction is represented over a finite time interval 
as a finite number of terms E01 



Then the frequencies and complex amplitudes are computed via an iterative 
procedure. For example, the first frequency is determined by the maximum 
amplitude of <p(a) 



where %(t) is an even-weight function. 

Numerical integration with a time step of 500 years over the time span of 
about 200 million years reveals that the inner planets of the Solar system are 
chaotic, due to the presence of two secular resonances, one due to Mars and 
Earth at 9 = 2{g^ — g%) — (S4 — S3) and another due to Mercury, Venus and 
Jupiter at a = (gi — g§) — (si — S2), where gi and Sj are the frequencies of the 
perihelions and nodes, respectively. 

Laskar's calculations revealed the chaotic behavior of Mercury's orbit with 
eccentricity variations up to 0.05, which results in its overlapping with the 
orbit of Venus and with inevitable escape of Mercury from its orbit. Chaotic 
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N 




(27) 




(28) 







50 100 
So , degrees 



Figure 1. The maximum, mean and minimum obliquity variations of the Earth depending 
on the initial obliquity eg. Note the chaotic zone at eg = 60° — 90° and the stability for 
its other values. The present obliquity of the Earth is within the stability zone due to 
the presence of the Moon, and transfers into the chaotic zone at the absence of the Moon 
(Laskar et al 1993). 



behavior was discovered also for the obliquities of the planets, particularly for 
Mars, varying from to 60 degrees. Obviously, this fact has to be taken into 
account while studying the past evolution of the atmosphere of Mars. Chaotic 
behavior of the obliquity of the Earth would range even wider, within and 85 
degrees, with all dramatic consequences for the climate of the Earth, however, 
only at the absence of the Moon. The Moon therefore, is dampin g th e obliquity 
variations up to 1.3 degrees, thus stabilizing the Earth's climate^. 

We see that, the chaotic effects are not only able to influence essentially 
the dynamics of the Solar system but even the Earth's climate. 
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Figure 2. The variation of the eccentricities of Mercury, Venus and Earth by time. The 
chaotic variations are suppressed for the Earth and Venus, but are significant enough for 
Mercury, and will lead to its escape from the Solar system within a period less than 3.5 
Gyr (Laskar 1994). 



4 Galactic Dynamics 

4-1 N-body gravitating systems and geodesic flows 

Many properties of statistical mechanics of globular clusters and galaxies can 
be studied considering a iV-body gravitating system described by Hamiltonian 



We will use a well known method existing in classical mechanics, the Mauper- 
tuis principlJ^, enabling one to represent a Hamiltonian system as a geodesic 




(29) 



(30) 



r a - n- 



(31) 
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Figure 3. The overlapping of the major semi-axes of orbits of Mercury and Venus due to 
chaotic effects (Laskar 1994). 



flow on some Riemannian manifold. In physical p roblem s this approach was 
firstly used by Krylo\El and in stellar dynamics ir F^ I ^ I 

By means of the Maupertuis principle the Hamiltonian equations 



dH dH 



dp^ 8r^ (32) 

are reduced to the geodesic equation 

V„u = (33) 
on the region of configurational space 

M = {W = E — V(r) > 0} 
with the Riemannian metric 

N 3 

ds 2 = [E- V(n, u rjv,3)] J2 E^) 2 ' ( 34 ) 

a—1 2—1 
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where E is the total energy of the system. The condition of conservation of 
the total energy of the system 

H{p,q) = E (35) 

is equivalent to the condition on the velocity associated with the geodesic 

|| u ||=1, (36) 

while the affinc parameter along the geodesic is determined by 

ds = V2(E - V(r))dt. (37) 

The statistical properties of the geodesic flow are determined from the 
Jacobi equation 

V„V„ti + Riem(n, u)u = 0. (38) 
For a vector field satisfying the orthogonality condition 

<n,u>=0, (39) 
the Jacobi equation can be written in the form 

rf2 \™ " 2 = -2K u , n || n || 2 +2 || V„n || 2 , (40) 

where K Uin is the two-dimensional curvature 

< Riem(n,u)u,n > < Riem(n,u)u 7 n > 

II u II II 71 II < u i n > II n II 

Jacobi equation has the solution 

II n(s) \\> \ || n(0) || cxp(V^2fc s ), s > (42) 

if 

fc = max{if u .„} < 0. (43) 

Then, the geodesic flow is an Anosov system. 

Thus the negativity of the two-dimensional curvature the criterion of the 
instability of the geodesic flow. 
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4-2 Chaos of spherical systems 

So, we have to calculate the two-dimensional curvature K Ut1l for the Hamilto- 
nia n ( ) fo r the gravitating N-body system. The Riemann curvature has the 

Riem^xup = --^[g^Wxp + gxpW^ - g pp Wxu - g^pW^x] 

[(g„xW v - gp V Wx)W p + {g yp Wx ~ gxpW^W^} 



where 



AW 2 

'^^[g^gxp - gp.xg»p\ II 9W 



gpv = WS^, 

dw d 2 W 



Q r pP j 3r^ dr u 

The analysis shows that K UtTl is sign-indefinite, which means that no universal 
function 

r : (all TV-body systems) — > R + 

exists and hence no unique relaxation time scale can exist for all TV-body 
gravitating systems. 

For spherical systems, however, the two-dimensional curvature at large 
N limit is shown to be strongly negative, since is determined by the scalar 
curvature R 

K u n = — - (44) 

' 3N(3N - 1) y ' 

which is negative for N > 2 

R - »i)(J- 5 L) ( v^- ( 3.- 1 ,-<o. <„, 

Based on the mixing properties of dynamical systems described in pre- 
vious sections, one can define the relaxation time for spherical systems. The 
explicit formula for that relaxation time taking into accou nt the n online ar in - 
teraction of all N bodies of the system has been estimated iJSEHl ( se e alsJ^I) . 
For real stellar systems its value is shorter than the two-body relaxation time 
scale, but longer than the dynamical (crossing) time. 

The disk gravitating systems can be studied using the Lie algebra of all 
vector fields with zero divergence on the two-dimensional torus T^" 1 ' , since 
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the kinetic energy of the element of the moving fluid induces a right-invariant 
Riemannian metric on SDif f(T 2 ). The principle of least action, which de- 
termines the motion of an incompressible fluid in terms of the geodesies of 
this metric, plays the role of the Maupertuis principle. One can show that, 
though the motion in disk galaxies is exponentially instable, the velocity field 
remains constant, so o ne c annot speak about a relaxation in the same sense 
as for spherical systemic 



4-3 Relative chaos is stellar systems 

The approach described above enables to consider instability of various con- 
figurations of stellar systems. Average the Jacobi equation over the geodesic 
deviation vector 



where 



zn, || h \\ 2 = f, 



and r u (s) denotes the Ricci curvature in the direction of the velocity of the 
geodesic (Ric is the Ricci tensor) 

Ric(u u) 3N ~ 1 
r u (s) = y — = ^ Kn^ujs), (n tl ±n„,n tl ±u). (47) 

The Ricci tensor has the expression 

Ricx P = ~ [AW 9V + (3N - 2)W Xp ] 

+ - " " "I dW "I 2 ■ 

Then the criterion of relative instability i 

the more unstable of two systems is the one with smaller negative r 

31\ 0<s<s, 

within a given interval < s < s„, i.e., this system should be unstable with a 
higher probability in the same interval. 

Numerical exploitation of the Ricci criterion of relative instability for the 
different models of stellar systems has shown that, e.g., a spherical system with 
a central mass is more unstable than a homogeneous one, spherical systems 
are more instable than disk-like ones, etcE^E2EE2 
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5 Galaxy Clusters: Substructure and Bulk Flows 



N-body gravitational systems, as we saw above, are exponentially instable 
systems. This fact gives a key to the possibility of reconstruction of certain 
properties of based on the limited observational information, which usually 
includes the 2D coordinates and ID (line-of-sight) velocities and the magni- 
tudes of the galaxies. We will show how one can reconstruct the hierarchical 
substru cture and the bulk regular flows of the subgroups in the clusters of 
galaxieJ3^l 

The developed S-tree technique is based on the geometrical methods of 
theory of dynamical systems discussed in the previous sections, namely, on 
the introduction of the concept degree of boundness of N particles. 

Consider two particles, so that x\{t) and x^if), t E (— T, T) are their 
trajectories when their interaction is taken into account, and yi{t) and y%{t), 
when the interaction is "switched off" . 

It is easy to see that the deviation of trajectories within certain time 
interval 

m = max M{xi(-)-yi(-)), (49) 

i— 1.2 

can be taken as a measure of the degree of boundness with respect to a local 
norm 

M{x{-)) = sup {\x(t)\,\x(t)\}. (50) 

t£(-T,T) 

Consider balls of radius r at each point of trajectories of the two inter- 
acting particles Xi . The union of those balls 

d(r)= |J B Bl(t) (r), i = l,2. 
te(-T.T) 

of such minimal radius m which contains all trajectories of the particles will 
denote the free corresponding particles. 

Two particles are considered to be p-bound for p > if m < p. 

This is easily generalized to any finite number of particles. N particles 
labeled by the set of integers A = {1, . . . , N} form a p-bound cluster if the 
distance between the corresponding trajectories of the system of interacting 
particles and free ones is less than the maximal deviation of all of the particles: 

m = maxA/"(x Q (-) - y a (-)) < p . (51) 
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One can then define the boundness function so that for the given local 
norm M 

V z (Y) = maxM(z a (-)-y a (-)), (52) 

where z a (t),y a {t) are the solutions of the systems of equations I and II re- 
spectively for some time interval (— T, T). In other words the boundness of 
y in Z is the maximum deviation of the trajectories of its particles taking 
into account only internal interactions compared to the situation when in- 
teractions with particles in Z are also included. Our goal is to split A into 
p-subsystems, i.e., to obtain the map £ for this choice of boundness function. 

The definition of a p-bound cluster given above can be reformulated now 
as a set of corresponding particles A being a connected subgraph of the graph 
r (equivalent the matrix) so that there is no other connected subgraph B 
including A: A C B. If one defines P as follows 

V X Y = max {D yz } , (53) 

z£X\Y 

the problem of the search of a p-bound cluster is reduced to that of a connected 
T-graph. 

The algorithm of the construction of tree-diagram based on the estima- 
tion of the two-dimensional curvature as containing information both on the 
coordinates and velocities of the all particles, is developed and applied to 
various clusters of galaxies. As a result subgroups of galaxies, "galaxy asso- 
ciations^, of specific dynamical properties, are detected in the studied Abell 
clusters^, triggering later observations by the provided lists of galaxies (see, 
e _gESj). The S-tree method, together with certain general assumptions on the 
velocity distribution function of galaxies can be used for the determination of 
the bulk velocities of the subgroups. 

Let us stress again, that these methods arise due to the nonlinearity of 
the N-body gravitating systems. 

6 General Relativity and Cosmology 

The impossibility of the direct application of results of theory of dynamical 
systems developed for Riemannian spaces is the main difficulty arising while 
studying the chaos in General Relativity and cosmology where one deals with 
pseudo-Riemannian spaces. Therefore, first, one has to reformulate the con- 
cepts described in previous sections and applied for astrophysical Newtonian 
systems, for the case of pseudo-Riemannian spaces. We will give the refor- 
mulation of the property of hyperbolicity and the covariant definition of the 
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Lyapunov exponents given i which are basic concepts for the study of 
chaos, and then, using it, will consider the stability of cosmological solutions, 
particularly, of inflationary ones. 

I will not discuss the mixmaster models which had essentially provoked 
the studies on chaos in cosmology, since they are covered in Kirillov's lectures 
at VIII Brazilian School of Cosmology and Gravitational For the further 
progress on those models in the context of Non-A belian ga uge, string theories 
and pre-Big-Bang scenarios I will refer to reviewJ^EIIEI 

6.1 Hyperbolicity in pseudo-Riemannian spaces 

Consider a geodesic flow on M, i.e. a group of mappings {<?*} of a space T X M 

T X M = {(x,u);x€ T x M,g(u,u) = \\u\\ 2 = A}, A = 0,±1. 

Each mapping performs a shift of a linear element £ = (x, u) along the geodesic 
on distance t. 

Let j(t) be a geodesic on M passing by a point x S M, and {E a } is a 
fixed n-dimensional basis on T X M. 

Transferring {E a } parallel along 7(t), i.e. getting a basis at every t, one 
has a Fermi basis on T X M . 

Each vector X S T 7 M can be represented via Fermi basis 



with the E-norm 



X(t) =X a (t)E a . 



\X\\ 2 E = X(Xa) 2 



for basis {E a }. 

Let {E a /} be another basis. Then a non-singular matrix <I>^ exists, such 

that 

Since both {E a } and {E a /} are Fermi bases, the latter relation has to be 
satisfied also for constant <I>^ . 
Then 

X b '(t) = Y 1 ^X a {t). 
In view of non-singularity of $^ , we can write 

CZ{X a ) 2 < ^{X b ,f < C- X E(X ) 2 
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or 

c\\x\\%<\\x\\%<c-'\\x\\l, 

where C is a positive constant. 

Definition of hyperbolicity. Geodesic j x (t) = S t (£ > ),\\j x (t)\\ 2 = A is A- 
hyperbolic, if there exist subspaces VK S (£*(£)) and W u (S t (^)) and W°(S t (^)) 
of the tangent space Tgt^T x M and numbers ^40,0 < /j, < 1, such that 

T st(c) T A M = VK s (5 t (0)^ u (5 t (e))VK (,5**(0), 

dS* r VT (£*(£)) = W s (S t+T (£)),dS T W u (S t (0) = W u (S t (0), 

where VF°(5'*(^)) is a ID space defined by the flow vector. 
For each t, r > and for a certain basis {E a } we have 

\\dSTv\\%<A*iP\\v\\%,veW'(S t ®), 
\\dSrv\\%>A-*»-*\\v\\%,veW«(S t {t)), 

where 

I Mil = WdwWl + \\Kv\\%,v G TT x M,ir x :TT X M^TM, 

and if is the mapping of connection V. 
The definition is {E a } invariant. 

Definition Geodesic flow is A— hyperbolic if its all geodesies are A- 
hyperbolic. 

Jacobi field is defined along the geodesic determined by the Jacobi equa- 
tion 

V u V u Y + R(u,v)Y = 0. 

Correspond now to each vector v e TT X M a solution Y(t) of Jacobi equation 
with initial conditions 

y„(0) = dnv,V u Y v (0) = kv. 

The resulting mapping 

f:v^Y v (t) 

is an isomorphism and 

dndS^v) = YvifyKd&v = V v Y v {t). 
From the Jacobi equation we have 

||^||i= \\Y v {t)\\ 2 E + \\W u Y u {t)\\ 2 E . 
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the latter equation and the Jacobi one enable to check the hyperbolicity con- 
dition. 

Definition. Lyapunov characteristic exponent for maximal geodesic 7 and 
vector v is defined as 

ln\\dStv\\% 
X(7, v) = hm sup -f . 

Definition. Geodesic 7 is stable if for any e > 0,3<5(e > such that 
from < 5 follows the condition ||cZ5^|| 2 < e for any t. Otherwise 7 is 

unstable. The latter two definitions are also basis-invariant. 

Let us now define a convenient basis. 

For arbitrary geodesic ^(t) we choose the following orthonormal basis at 
point 7(t) 

Ea = 7(0) = u,E u ...,£„_i; 
g{E a ,E b ) = 

where E a is a dual basis. 

If the following conditions are satisfied 

V u E a = = V U E\ 

then the basis on T°M can be defined as 

E a = u,g(E a ,E b ) = 



and on T X M 



For the vector field 



E = u,g(E a ,E b ) 



Y(t) = Z a (t)E a 

the Jacobi equation can be written in the form 

Z a (t) + KZ(t)Z b (t) = 0, 

where 

Kl =< E a , R(u, E~b)u >= R a bcd u c u d . 
For the above defined basis on T°M and the Jacobi equation we have 

Z = 1. 

which means that none of geodesic flows can be 0-hyperbolic. 
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It can be shown that the definitions given above for spaces of Lorentzian 
signature (-,+,...,+), can be generalized for the signatures (-,,-, +,,+). 

The covariant definition of hypcrbolicity and Lyapunov exponents given 
above enable the consideration of the stability problem of cosmological solu- 
tions. 



6.2 The ADM principle and geodesic flows in Wheeler- DeWitt 
superspace 

The problem of the stability of cosmological solutions is a problem of sta- 
bility in Wheeler-DeWitt superspace. We will consider this problem using 
the method of geodesic flowg32l. So, first, we have to define the Hamilto- 
nian system, then reduce it to a flow of geodesies using the definition of the 
hyperbolicity given in the previous section. 

Arnowitt-Deser-Misner (ADM) method provides the scheme of the sought 
Hamiltonian formulation, assuming as given the 3-geometries of the initial and 
final Cauchy hypersurfaces. We will consider locally isotropic and homoge- 
neous cosmological models with scalar field when the metric can be given 
as 

where 

a 2 = ^_ [ J - h i/2 d z x] -\ h = defh l3 
We consider the Lagrangian of the scalar field 

and the action 




p a da + p x dx - ATUDMdt 



where the ADM Hamiltonian is 

Hadm = ^[-p* + p£] + e 3 «[U(x) - ±e~ 2 % 

where 

4 G" ^ 

a = Ina, \ = o~(j), U(x) — V{4>). 

o 

As usual variation with respect the lapse function N leads to the condition 

Hadm = 0. 
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To reduce the Hamiltonian system to the geodesic flow let us split the hyper- 
surface into the following regions 

W + {x\V(x) > 0}, W~{x\V(x) < 0}, 
so that if the metric in region W~ is Riemannian, then 

9a b v a v b = -2V > 0, dr = {g ab u a u b / - 2V)~^ 2 ds 
and we can write the variation 

extI\ H=0 = ext j {-2V) 1/2 {g ab u a u b ) 1 ' 2 ds = ext J {G ab u a u b ) 1 ' 2 ds. 

Here G ab = —Vg ab is also a Riemannian metric. 

Choosing the affinc parameter s in order to satisfy the condition 

\\u\\ 2 = G ab u a u b = 1 

we have 

G ab u a u b = -Vg ab v a v b (dT/ds) 2 = 2V 2 (dT/ds) 2 . 

Rcparameterizing the afline parameter 

ds = 2^ 2 (-V)dT, 

we arrive at the flow of geodesies in the region W~ 

H = l/2g ab Pa Pb + V - {G ab = -Vg ab , ds - 2 1 / 2 (-V)dr, ||u|| 2 = 1}. 

As regards for the region W + , the classical system cannot end up in it. 

Now, if the metric g is pseudo-Ricmannian, following the same scheme we 
end up with the geodesic flow 

H^{|V|g ab , 2 1 / 2 |V|dr,-signV}. 

Thus, we reduced the ADM Hamiltonian system to a geodesic flow on a 
pseudo-Ricmannian manifold. 

To study the stability of the geodesic flow we have to proceed from the 
Jacobi equation, which has the form 

d 2 dz^ 
^ +7 (z)_+ „}(*) = (), 

where 

7 = -^ln\V\,u) = 2V 2 K},i,j,= l,2...k - 1. 
Using the variables 

z l = AY 1 
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we arrive to the equation 

Y l + 2(A/A + 7 )r* + [{A/A + iA/A))8) + w^Y* = 

where 

A/A = -l/2j,A = \V\ 1/2 . 
Particularly for the case k = we have the simplified expressions 

G = e 6a \U\,\V\ = e 3a \U\, 

. . da 1 dU 

^ T) ^- 3 d^-u^ 

, N 1 DlnG UnD\V\ 
K{t) 



2 G 2 e 6a |[/| ' 

a 2 d 2 

uj{T ) = UDln\U\,D = -^ + — 2 . 

6.3 Stability of inflationary solutions 
Consider the scalar field 

U = X X n /n 

and the conditions 

H«H 2 ,± 2 «U,\x/x\ «H 
and x >> 1 at a — > oo. Then we have 

and the Jacobi equation 

^-[-^ + (|ff) 2 ]^ = o. 

A Z 

From the Einstein equation 

H + 'iH 2 = 6U, 

and in view of the condition H « H 2 we have H 2 ~ 2t7\ If x >> (2n/3) 1/2 , 
then nU/x 2 « 9/4_ff 2 , and we have the Jacobi equation in a simple form 

Y- (9/4a 2 )Y~ = 0. 
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From its solution we obtain 

u' 4 

S = z +z c± const U + const -jjj ~ const x™ + const x 2n 4 - 

We see that S decreases for any n > 2 and therefore we have Lyapunov stability 
of the inflationary solutions. The last formula enables also to obtain the law 
of the decay of perturbations at various n. For example at n = 4 we have 
exponential decay of perturbations, and the larger is n, the more stable the 
solution is. 

Thus we showed how one can deal with the stability problem in pseudo- 
Riemannian spaces, and illustrated this on inflationary solutions. 



7 Instability in Superspace 

How typical is the given cosmological solution? This is a basic question posed 
since the early days of the study of the Einstein equations. In the context 
of later developments, particularly in quantum cosmology, the question can 
be reformulated in the form: to what degree are the minisuperspace models 
typical in superspace given the huge extrapolations involved? The considera- 
tion of perturbed minisuperspace models by Hawking and other authors still 
involves extrapolation in the absence of a deeper merging of quantum theory 
and gravity. 

The study of dynamics in the Wheeler-DeWitt superspace, more precisely, 
the properties of geodesic flows in superspace can provide a more general view 
of how typical the minisuperspace models are. The problem, however, is far 
more difficult than the one posed in conventional hyperbolicity theory, since 
one deals both with infinite dimensional and pseudo-Riemannian manifolds. 

However, hyperbolicity can be defined for such man ifolds and we will 
consider the case of homogeneous cosmological models ™l 

The metric of Wheeler-DeWitt superspace is 

Qijkl _ i[ykyJ _|_ ryil^jk _ 2-f i 3j kl ] l 

where 

7 = defrfij,i,j, k,l = l, ...,n. 

Then one can derive 
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where 

£ = 4((n - l)/n) 1 /2 7 V4 ,,4,5 = 1, ...,n(n + l)/2 - 1 , 
Then, moving to the subspace of superspace 

^ = ()ii;7ii^,7 = i) 

with a metric induced by the metric of the superspace W 

G AB dd A de = *r(7- 1 07/0£V 1 fry/0£ B ) de A * B > 

the existence of non-zero Lyapunov numbers 

YsiUJi = -n/A < 
can be shown for the solutions of the Jacobi equation 

a? exp[(±(- Wi ) 1/2 t]. 

This implies the exponential instability of the geodesic flow in that subspace 
of the superspace. 

For models with a scalar field Armen Kocharyan^l was able to show that 
the instability is exponential if: 

1. Gravitational and matter fields vary quickly with respect the potential; 

2. The Universe undergoes inflation in a local domain. 

The smaller arc the dimension and the number of scalar fields, the stronger 
is the instability. 

These results lead to the following general conclusion: 

The quantized system in a finite-dimensional submanifold is not 
typical to that in superspace due to the existence of virtual pertur- 
bations along the frozen directions which are unstable. 

This implies that minisuperspace models cannot be considered as fair 
approximations of superspace models. 

8 Conclusion 

I stop here. As we saw, chaos is an inevitable ingredient of the Universe, but 
it needs particular efforts to deal with. 

The spectrum of the methods applied will obviously grow further. 1 con- 
clude with mentioning the use of Kolmogorov complexity (algorithmic in- 
formation) for the study of properties of Cosmic Background Radia tion and 
relations of the latter with thermodynamical and cosmological 
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